Analysis of Kelly-optimal portfolios 
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, Abstract 



We investigate the use of Kelly's strategy in the construction of an optimal portfolio 
of assets. For lognormally distributed asset returns, we derive approximate analytical 
results for the optimal investment fractions in various settings. We show that when mean 
CO ' returns and volatilities of the assets are small and there is no risk-free asset, the Kelly- 

optimal portfolio lies on Markowitz Efficient Frontier. Since in the investigated case the 
Kelly approach forbids short positions and borrowing, often only a small fraction of the 
available assets is included in the Kelly-optimal portfolio. This phenomenon, that we call 
condensation, is studied analytically in various model scenarios. 



1 Introduction 



The construction of an efficient portfolio aims at maximising the investor's capital, or its 
return, while minimising the risk of unfavourable events. This problem has been pioneered by 
^ ■ Markowitz in [l|, where the Mean- Variance (M-V) efficient portfolio has been introduced: it 

minimizes the portfolio variance for any fixed value of its expected return. Since this rule can 
be only justified under somewhat unrealistic assumptions (namely a quadratic utility function 
^Sj ■ or a normal distribution of returns, in addition to risk aversion), it should be considered as a 

, first approximation of the optimisation process. Later, several optimisation schemes inspired 

by Markowitz's work have been proposed 0, 0, 0]- For a recent thorough overview of the 
portfolio theory see |5|. 

A different perspective has been put forward by Kelly in where he shows that the 
optimal strategy for the long run can be found by maximising the expected value of the 
^ ' logarithm of the wealth after one time step. The optimality of this strategy has long been 

treated and proven in many different ways 3, S, 0], according to [10], it was successfully used 



in real financial markets. For an overview of its continuous time limit see [11[. Recently, 
the superiority of typical outcomes to average values has been discussed from a different 
point of view in [l^, [13] ■ Although the Kelly criterion does not employ a utility function, as 
pointed out by the author himself, a number of economists have adopted the point of view of 
utility theory to evaluate it [ij, |l3, [3, [itIi. Various modifications, such as fractional Kelly 
strategies [18| and controlled drawdowns [l9| , have been proposed to increase security of the 
resulting portfolios. A thorough review of the advantages, drawbacks and modifications of 



the Kelly criterion is presented in 20(]. For an exposition of the Kelly approach in the context 



of information theory see 21 ] . 



In this paper, we shall discuss the original Kelly strategy in the framework of a simple 
stochastic model and without assuming the existence of a utility function. We will present 
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approximate analytical results for optimal portfolios in various situations, as well as numeri- 
cal solutions and computer simulations. We will show that, in the limit of small returns and 
volatilities, when there is no risk- free asset, the Kelly-optimal portfolio lies on the Efficient 
Frontier. Furthermore, we shall analytically study the conditions under which diversification 
is no longer profitable and the optimal portfolio "condensates" on a few assets. Such con- 
densation (or underdiversification) is said to be typical for the Kelly portfolio [2(1] and here 
we examine it in various model scenarios. Finally, we will consider the fluctuations of the 
logarithm of wealth as a measure of risk, and compare it with the classic M-V picture. 

This paper is organised as follows. After introducing a multiplicative stochastic model for 
the dynamics of assets' prices, we briefly list the main results of the Markowitz Mean- Variance 
approach. In Sec. [21 we apply Kelly's method to our model, analysing the case of one, two 
and many risky assets, both with and without additional constraints. Finally, a combination 
of the Markowitz Efficient Frontier with the Kelly strategy is investigated. In the appendix 
we explain the approximations used in this paper as well as a generalisation of the model to 
the case of correlated asset prices. 



1.1 A simple model 

We shall study the portfolio optimisation on a very simple model which leads to lognormally 
distributed returns. Consider assets, whose prices (i = 1, . . . , N) undergo uncorrelated 
multiplicative random walks 

Here the random numbers r]i{t) are drawn from Gaussian density distributions of fixed mean 
rrii and variance Di, and are independent of their value at previous time steps. This model 
can be easily generalised to the case of non-Gaussian densities and correlated price variations 
as it is discussed in appendix iBj the influence of correlations on the Kelly portfolio is investi- 
gated in [i^l . We assume that the investor knows exact values of the parameters rrij , Di — for 
the effects of wr ong pa rameter estimates and the details of the Bayesian parameter- learning 



20 . l23l . |2j] . We further assume the existence of a risk- free asset paying zero 



process see e.g. 
interest rate. 

For the sake of simplicity, we do not include dividends, transaction costs and taxes in 
the model. Hence, the return of asset i is Ri{t) := [pi{t) — pi{t — 1)] /pi{t — 1) = e''*^*^ — 1 is 
lognormally distributed with the average fii := E{Ri) = exp[mi -|- Di/2\ — 1 and the volatility 
af := E{{Ri — fJ-i)'^) = {exp[Di] — 1) exp[2mj -|- Di]. With E we denote averages over the noise 

A portfolio is determined by the fractions Qi of the total capital invested in each one of 
available assets; the rest is kept in the risk-free asset. Since and Di are fixed, both 
the Kelly strategy and the Efficient Frontier use one time step optimisation and the basic 
quantity is the wealth after one time step Wi. If we set the initial wealth to 1, Wi has the 
form 

N 

Wl = l + Y,Q^Ri = ^ + Rp, (2) 

i=l 

where Rp := ^i^iQiRi is the portfolio return. To simplify the computation we assume 
infinite divisibility of the investment. Thus, the investment fractions qi are real numbers and 
do not need to be rounded. 

In the portfolio optimisation, some common constraints are often imposed and can as well 
be applied in the present context. For instance, the non-negativity of the investment fractions 
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Figure 1: The expected return /ip versus the standard deviation of the portfoho ap. The 
assets are described by the following parameters: m\ = 0.1, Di = 0.04, m2 = 0.15, D2 = 0.09, 
7713 = 0.2, 1)3 = 0.25 (in the graph they are shown as full circles). The dashed line represents 
the CML from Eq. ([3]), the solid line is Efficient Frontier given by Eq. (I5|), the tangent point of 
the two is the Market Portfolio. The thick part of EF marks the region where all investment 
fractions are positive. 

> forbids short positions. The condition Yld=i Qi = ^ indicates the absence of a riskless 
asset and J2f=i Qi ^ ^ does not allow the investor to borrow money. 

1.2 The Mean- Variance approach 

The unconstrained maximisation of the expected capital gain results in the investment of 
the entire wealth on the asset with the highest expected return; this strategy is sometimes 
referred to as risk neutral. If the investor has a strong aversion to risk, on the other hand, 
one might be tempted to simply minimise the portfolio variance ap = Xli^o ^f'^f- This leads 
to invest the entire capital on the risk-free asset with no chance to benefit from asset price 
movements. The Mean- Variance (MV) approach is much more reasonable as it allows to 
compromise between the gain and the risk. Here we recount basic results of this standard 
tool. 

With the desired expected return fixed at E{Rp) = fip, the constrained minimisation of 
the portfolio variance ap is performed using the Lagrange function C = E{R],) +-i{E{Rp) - 
fip) with a Lagrange multiplier 7. The resulting optimal fractions are 



For Hp = 0, (ji = for all assets. As we increase fip, all optimal fractions qi grow in a uniform 
way and their ratios are preserved. At some value ^p we reach Y^iji = 1, which means we 
are investing the entire capital. Any further increase would require to borrow money, with 
Eq. ([3]) remaining valid as long as the borrowing rate equals the lending rate (both set to zero 
here). The relation between ap and /ip is 




(3) 




(4) 



This equation is often referred to as Capital Market Line (CML). 
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If there is no risk- free asset in the market, one has to introduce the additional constraint 
Si^i = 1- It fohows that 

2 _ Cpfip — 2Cifj,p + C2 , . 

— 777; 7^2 ■ 

The functional relation between the optimised ap and fip is called Efficient Frontier (EF). 
Since there is only one point on the CML where - qi = 1, this line is tangent to the EF. The 
results of this section are plotted in Fig. [T] for a particular choice of three available assets. 



2 The Kelly portfolio 

When the investor's capital follows a multiplicative process, after many time steps is its ex- 
pected value strongly influenced by rare events and in consequence it is not reasonable to 
form a portfolio by simply maximising E{W{t)). The Mean- Variance approach tries to solve 
this problem in a straightforward, yet criticisable way. We support here the idea that an 
efficient investment strategy can be found by maximising the investment growth rate in the 
long run, which is, under the assumption of fixed asset properties, equivalent to maximising 
the logarithm of the wealth Wi after one time step 0]. Thus the key quantity in the con- 
struction of a Kelly-optimal portfolio is v := E{\b.Wi)^ the average exponential growth rate 
of the wealth. We remind that the quantity In Wi is not a logarithmic utility function. 

In V is optimised in a similar context and the authors claim that their procedure 
corresponds to maximising the median of the distribution of returns. They consider short 
time intervals and thus small assets returns. Assuming Rp <^ 1 (very small portfolio return), 
they use the approximation ln(l + Rp) ~ Rp — -Rp/2 of the logarithm in the expression of 
V before maximising it. However, while such an expansion is only justified for Rp <C 1, the 
maximum of the resulting function is at Rp = 1, in contradiction with the hypothesis. We 
will develop a different approximation in the following. 

First, the unconstrained maximisation of v is achieved by solving the set of equations 
dv/dqi = (i = 1, . . . , N). After exchanging the order of the derivative and the average, we 
obtain the condition 

In our case, Ri has a lognormal distribution and to our knowledge, this set of equations cannot 
be solved analytically. With the help of the approximations introduced in the appendix we 
shall work out approximative solutions for some particular cases. We emphasize an important 
restriction which applies to all solutions of Eq. ([6]). Since returns Ri lie in the range (—1, 00), 
when X^i^i 9i > 1 or when there is an investment fraction qi < 0, there is a nonzero probability 
that Wi is negative and hence v = E{ln Wi) is not well defined. Since Kelly's approach focuses 
on the long run, it requires strictly zero probability of getting bankrupted in one turn. As 
a consequence, for lognormally distributed returns any Kelly strategy must obey % > and 
X^i^i < 1, i-e. both short selling and borrowing must be avoided. 



2.1 One risky asset 

Let us begin the reasoning with the case of one risky asset. We want to find the optimal 
investment fraction q of the available wealth. The remaining fraction 1 — q we keep in cash 
at the risk-free interest rate which, without loss of generality, is set to zero. This problem 
is described by Eq. Q in one dimension; even this simplest case has no analytical solution. 
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Figure 2: The optimal portfolio fraction q: a comparison of the analytical result Eq. ([7]) with 
a numerical maximisation of E(\x].Wi). 



Nevertheless, for a given D, one can ask what is the value m< for which it becomes profitable 
to invest a positive fraction of the investor's capital in the risky asset. This can be found 
imposing g = in Eq. (I6|), yielding m< = —D/2. Similarly, the value m> for which it becomes 
profitable to invest the entire capital can be found by imposing ^ = 1, yielding m> = D/2. 

We shall look for approximate solutions that are valid for small values of which is the 
case treated in appendix El Using approximation Eq. ([20]) in Eq. ^ gives 

+ ^ hr^ = 0. 



l-q + qe"^ 2 {I - q + qe^ f 

With respect to q, this is merely a quadratic equation. Since the solution is rather long, we 
first simplify the equation using m, D <C 1 as in Eq. (f2T]l . leading to the result 



1 "T- /r,N 

Since borrowing and short selling are forbidden, for m < —D/2 is g = and for m > D/2 
is g = 1. When asset prices undergo a multiplicative random walk with lognormal returns, 
both m and D scale linearly with the time scale and hence q does not depend on the length of 
the time step. Notice also that substituting m = ^D/2 gives g = and g = 1, in agreement 
with the bounds we found before by exact computation. The first order correction to Eq. ([7]) 
is m{Am? — D'^)/AD'^ which is, for m e [—D/2, D/2], of order 0{m). The validity of the 
presented approximations can be easily tested by a straightforward numerical maximisation 
of E(\a.W\). As can be seen in Fig. [21 the numerical results are well approximated by the 
analytical formula Eq. ^ even for D = \. 

Notice that, for m,D <^ 1, one can approximate ~ m + D/2 and cr^ ~ D, which 
makes the optimal portfolio fraction derived above equal to q' = + 1^^) obtained in [13]. 

However, if we check the accuracy of q' , we find a relative error up to 3% for D = 0.01, and 
for D = 0.25 we are already far out of the applicability range with an error around 50%. Also, 
Eq. ([?]) is for m, L> <C 1 identical to the classical Merton's result [iH] which is derived under 
the assumption of continuous-time, non-zero consumption of the wealth, and a logarithmic 



utility function. In 11, |26[], the result q = ^ + m/D is derived for the continuous time limit 



of our model: m, D — > 0, m/D = const. 
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2.2 Constrained optimisation 

The optimal portfolio fractions qi can be derived from Eq. ^ also for > 1. Using the same 
approximations as in the single asset case, we obtain the general formula qi = 1/2 + rrii/Di 
for i = 1, . . . , N . In our case, Kelly's approach forbids short selling and hence the assets with 
gi < do not enter the optimal portfolio. Since borrowing is also forbidden, if X^^^ qi > 1, 
we have to introduce the additional constraint qi = 1. This can be done by use of the 

Lagrange function 7) = i' + 7 ( X]j=i ~ ^) where q is the vector of investment fractions. 
The optimal portfolio is then the solution of the set of equations 

^ R 



where i?j = e''* — 1. Using the same approximations again, one obtains the general result 

1 +7 

The Lagrange multiplier 7 is fixed by the condition X^jLi 9i = 1- It can occur that even 
a profitable asset with rrii > —Di/2 has a negative optimal investment fraction. Since in 
our case the Kelly approach forbids short selling, this asset has to be eliminated from the 
optimisation process. In consequence, under some conditions, only a few assets are included 
in the resulting optimal portfolio. This phenomenon, which we call portfolio condensation, we 
study closer in sections 12. 31 and 1 2.41 An alternative approach to the constrained Kelly-optimal 



portfolio is provided by the Kuhn- Tucker equations (see 2l|) which, however, can be shown 
to be equivalent to Eq. ([8]). 

Now we can establish an important link to Markowitz's approach: in the limit /ii,(Tj — > 
the Kelly portfolio lies on the constrained Efficient Frontier (no short selling allowed). We 
shall prove this statement in the following. When all the assets have small /Xj and cTj , in Eq. ([3]) 
and Eq. ^ we can approximate fii ^ nii + Di/2 and af ^ Di, leading to the approximative 
relation for the Efficient Frontier 

(7o,^-2C,,p + C. where C. = fi!^^i±^. (10) 
^ C0C2-C2 Di ^ ^ 

For the Kelly portfolio we need to work out a similar approximation. Using the condition 
= 1, for 7 in Eq. dH) we obtain 7 = (1 — Ci)/Co. In the relations fip = YliQif^i ^^'^ 
ap = "^iqfcrf we use the approximations for fii,ai introduced above. After substituting qi 
from Eq. ([9]), for the Kelly optimal portfolio we get 

C0C2 — Cf + Ci 2 C0C2 — Cf + I 

IJ-K = , o-x = • (11) 

(-^0 L'O 



Both in Eq. (jlOp and Eq. (jlip we consider only the assets that have positive investment 
fractions. Now it is only a question of simple algebra to show that ^ik and ax given by 
Eq. (jlip fulfill Eq. (jlOp . which completes the proof. Similar, yet weaker, results can be found 
in the literature. For instance, Markowitz states in [itI] that "on the EF there is a point which 
approximately maximizes E{[nWi)y 

Obtained results are illustrated in Fig. [31 where we plot the Efficient Frontier, the con- 
strained Efficient Frontier, and the Kelly portfolio for the same three assets as in Fig. [1] 
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Figure 3: The Efficient Frontier (EF, thin sohd line), EF without short selhng (bold line), 
and the Kelly portfolio (solid square) in a particular case of three assets (asset parameters as 
in Fig.ll]). 

While the original EF is not bounded (for any fip exists appropriate crp), the constrained EF 
starts at the point corresponding to the full investment in the least profitable asset and ends 
at the point corresponding to the most profitable asset. The two lines coincide on a wide 
range of fip. In agreement with the previous paragraph, the Kelly portfolio lies close to the 
constrained EF. 



2.3 Condensation in the two asset case 

To illustrate the condensation phenomenon we focus on a simple case here: two risky assets 
plus a risk-free one, borrowing and short selling forbidden. As we have already seen, without 
constraints qi = 1/2 + rrii/Di. Therefore, when rrij < —Di/2, qi is negative and due to 
forbidden short selling, asset i drops out of the optimal portfolio. In Fig. H] this threshold is 
shown for i = 1, 2 by dashed lines. In the lower-left corner (A) we have the region where both 
assets are unprofitable and the optimal strategy prescribes a fully riskless investment. 

When the results of the unconstrained optimisation sum up to one {qi -|- ^2 = 1)) we are 
advised to invest all our wealth in the risky assets. If both assets are profitable, this occurs 
when mi/ Di+m2/ D2 = 0. When only asset i is profitable, we should invest the entire capital 
on it only when rrii equals at least Di/2. In Fig. H] these results are shown as a thick solid 
curve. 

Since borrowing is not allowed, in the region above the solid line constrained optimisation 
has to be used. The condensation to one of the two assets arises when the optimal fractions 
((71,(72) are either (1,0) or (0, 1); we can find the values m'^^ and m'/ when this happens. By 
eliminating 7 from Eq. ([8]) and substituting (?i = 1 and 52 = we obtain the condition for the 
condensation on asset 1: E[{eP^ — l)/e''^) = £'((e^^ — l)/e''^). This can be solved analytically, 
yielding 

, D1 + D2 



mi 



7712 



(12) 



This equation holds with interchanged indices for the condensation on asset 2, thus m'/ = 
ra2 — {Di-\-D2) /2. Finally, for m'l < mi < m\ the optimal portfolio contains both assets. The 
crossover values m'^ and m'2 are shown in Fig. [Has dotted lines. They delimit the region where 
the portfolio condensates to only one of two profitable assets. A complete "phase diagram" of 
the optimal investment in the two assets case is presented in Fig. H] for a particular choice of 
the assets' variances. Interestingly, growth-rate optimising strategies have their importance 
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Figure 4: The phase diagram of the two-asset system with Di =0.1 and D2 = 0.2. In region 
A the investor is advised to use only the risk- free asset. In regions Bi, B2, and C the optimal 
investment is still partially in the risk- free asset. In regions Di, D2, E, Fi, and F2 one should 
invest everything in the risky assets. While in regions C and E the investment is divided 
between the two assets, in shaded regions Fi and F2 a nontrivial condensation arises: one is 
advised to invest all wealth in one asset although the other one is also profitable. 



also in evolutionary biology [27] where a similar condensation phenomenon has been observed 
when studying evolution in an uncertain environment (2^ . 

2.4 Many assets with equal volatility 

We investigate here the case of an arbitrary large number N of available assets where Kelly's 
approach, which forbids borrowing and short selling, gives rise to a portfolio condensation. 
While the optimal portfolio fractions are given by Eq. to find which assets are included in 
the optimal portfolio is a hard combinatorial task. To obtain analytical results, we simplify 
the problem by assuming that the variances of all assets are equal, Di = D (i = 1, . . . , N). 
The number of assets contained in the optimal portfolio is labeled as M and the assets are 
sorted in order of decreasing rrii {mi > m2 > ■ ■ ■ > mi\j). 

If the unconstrained optimisation does not violate forbidden borrowing, then profitability 
of an asset (i.e., rrii > —D/2) is the only criterion for including it in the optimal portfolio. 
When the constrained optimisation is necessary, the optimal portfolio is formed by starting 
from the most profitable asset mi , and adding the others one by one until the last added asset 
has a nonpositive optimal fraction qm+i ^ 0. Summing Eq. ([9]) from 1 to M, we can write 7 
as 7(M) = D{M~^ — 1/2) — X^i^ -^o'^ ^ given realisation of {mi}, we can obtain the 
resulting portfolio size by finding the largest M that satisfies qm > 0, which leads to 



This relation tells us how many assets we should invest on, once their expected growths and 
volatility are known. Notice that for M = 2 and Di = D2 = D, this result is consistent with 
that of Eq. (I12p where a special case of the condensation on two assets is described. 

Let us follow now a statistical approach. If all are drawn from a given distribution 
/(m), the value of M depends on the current realisation. The characteristic behaviour of 





(13) 



8 



the system can be found by taking the average over all possible realizations and replacing nii 
by rrii. The resulting typical portfolio size Mt captures this behaviour and depends on the 
distribution f{m) and on the number of available assets N. 



2.4.1 Uniform distribution of m 

Let us first analyse the case of a uniform distribution of mj within the range [a, b]. First we 
assume that all assets are profitable, i.e. a + D/2 > 0. For rrii are sorted in decreasing order, 
one can show ffii = b — {b — a)i/{N + 1). Since mj declines with i linearly, according to Eq. ^ 
so does Qi. Substituting rrTj for rrij in Eq. (fT3l) and replacing > with =, we can estimate the 
typical number of assets in the optimal portfolio Mt- Assuming Afj- ^ 1, the solution has 
the simple form Mt = yj2ND/{b - a). 

We are now able to generalise this result to the case where not all assets are profitable, i.e. 
rrii + D /2 < Q for some Vs. In the extreme case b + D/2 < and all assets are unprofitable, 
leading to M = Mt = 0. The opposite extreme is realized for mM^+i + -D/2 > which falls 
in the previously treated case because the number of profitable assets is larger than Mt- In 
the intermediate region, only the assets with nii > —D/2 are profitable and enter the optimal 
portfolio. On average, they are N{b + D/2)/{b — a). All together we have the formula 



Mt 



(6 + L>/2<0), 

{b + D/2>0, fnMr + D/2<0), (14) 



^ {mM, + D/2>0). 



In Fig. O left, an illustration of a particular system (a = x — L, b = x + L, x = —0.05, 
motivated by [2^) is shown. We plot Mt and /ip as functions of L. When L is small, 
all available assets are unprofitable and the optimal strategy is to keep the entire capital 
at the risk-free rate. As soon as the first profitable assets are added to the system, the 



optimal portfolio includes all of them, until it saturates at the value \JND jL (in Eq. ([T] 
we substitute b — a = 2L) . A further increase of L widens the distribution of m and enlarges 
the gaps between profitable assets. It becomes, as a consequence, more rewarding to drop the 
worse ones and Mt decreases. The analytical solution, displayed in Fig. [5] as a solid line, is 
in a good agreement with the numerical results (shown as symbols). Although no single-asset 
portfolio arises in this case, the relative portfolio size is Mt/N ~ l/-v/iV and hence in the 
large limit, the optimal portfolio condensates to a small fraction of all available assets. 

A more flexible measure of the level of condensation is the inverse participation ratio, 
defined as 7^ = l/X^i^i^f- It estimates the effective number of assets in the portfolio: 
when all investment fractions are equal, TZ = N, while when one asset covers 99% of the 
portfolio, TZ ^ 1. Concerning the typical case, using Eq. ([9]) we can write = A — Bi, 
B = {b — a)/[D{N + 1)], the detailed form of A is not needed for the solution. We assume 
that the number of profitable assets is larger than the typical size of the optimal portfolio 
Mt- Consequently, passing from i = 1 to i = Mt, Qi decreases linearly to zero and we can 
use the identity E,=i(^ - Bif = EZliBif to obtain 

7^ = [B^Mt{Mt + l)(2Mr + l)/6]"^ « [b'^M^/S]'^ « fMy. (15) 

In the last step we used Eq. (I14p for the typical size of the condensed optimal portfolio. We 
see that the uniform distribution of mj leads to the inverse participation ratio proportional 
to the number of assets in the portfolio. In the right graph of Fig. [SI Eq. p5|) is shown to 
match the numerical solution (based on Eq. ()13p ) for various numbers of available assets. 
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Figure 5: Left: In the case with rrii uniformly distributed in [x — L,x + L] we plot the size 
of the optimal portfolio Mt (solid line represents analytical result Eq. (jl4p . symbols are 
numerical results averaged over 10 000 repetitions). With the dashed line the portfolio return 
fip (measured in percents) is shown. The parameters are = 1000, D = 0.01, x = —0.05. 
Right: In the case of N assets with D = 0.01 and mj uniformly distributed in [0, 0.1], we plot 
the average size of the optimal portfolio Mt, and the inverse participation ratio TZ. Solid 
lines are the analytical solutions reported in Eqs. (fT^ and (fT3]) . symbols stand for numerical 
results. 

2.4.2 Power-law distribution of m 

Now we treat the case of a distribution /(m) that has a power-law tail: f{m) = Cm~°'~^ for 
m > mmin- As long as M <C A^, the properties of the assets included in the optimal portfolio 
are driven by the tail of f{m). In consequence, the detailed form of f{m) for m < m^i^ is not 
important here. We assume that only a fraction r of all assets falls in the region m > mmin- 
Instead of seeking the typical portfolio size Mt, we shall limit ourselves to finding the 
conditions when a condensation on one asset arises. With this aim in mind, we put M = 2 in 
Eq. (jl3p . obtaining the equation mi — m2 = D. When mi — m2 > D, only asset 1 is included in 
the optimal portfolio. By replacing mi and m2 with their medians mi and m2, one obtains an 
approximate condition for a system where such a condensation typically exists. Using order 
statistics [s^ we find the following expressions for the medians: mi = mmin (-^^^Z In 2)^/", 
'7^2 = mmin(A^'"/l-68)^/". The equation mi — m2 = D thus achieved can be solved numerically 
with respect to a. In this way we find the value ai below which the optimal portfolio 
typically contains only the most profitable asset. In Fig. [6] we plot the result as a function 
of D. For comparison, the outcomes from a purely numerical investigation of the equation 
P(mi — m2 > D) = 0.5 are also shown as filled circles. Our approximate condition has the 
same qualitative behaviour as the simulation, showing that the use of median gives us a good 
notion of the optimal portfolio behaviour. 

3 Efficient frontiers 

Markowitz's Efficient Frontier is the line where efficient portfolios are supposed to lie in the 
Mean- Variance picture. Here we would like to follow the same procedure, using typical instead 
of average quantities. To capture a typical case, we replace the portfolio return Rp by InW^i 
in all averages of Sec. 11.21 According to the formula E{{x — E{x))'^) = E{x'^) — E[x)'^ we can 
minimise E{{lnWi)'^) - E{lnWif instead of E{{lnWi - E{lnWi)f). With the constraints 
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Figure 6: Values of the power law exponent ai at which condensation to one asset arises for 
N = 1000, r = 0.1, and 

'^min — 0.1. We compare the semi-analytical result obtained by 
solving fhi — 1712 = D (shown as the solid line) with a numerical simulation of the system 
(shown as filled circles). Below the line, the optimal portfolio typically contains only one 
asset. 



E{lnWi) = vp and "^^^i qi = 1, the Lagrange function has the form 

AT 

£ = E{{\nWif)+^i{E{lnWi) - ^^p) + 72 ( 9* - l) • (16) 

i=l 

Its analytical maximisation leads to complicated equations and thus it is convenient to inves- 
tigate the system numerically; we do so in the particular case of three assets used in Fig. [1] 
Due to the two constraints there is effectively only one degree of freedom for the minimisation 
of E{(\nWi)'^) and the numerical procedure may be straightforward. 

For the resulting portfolios we can compute expected returns and variances which allows 
us to add this "Logarithmic Efficient Frontier" (LEF) to the ap-fip plane depicted in Fig. [TJ 
The result is shown in Fig. [71 where the solid line is again Markowitz's EF. Solid circles 
correspond to the three individual assets. The dashed line represents LEF (obtained by the 
numerical optimisation described above) and the thick gray curve is the region where both 
EF and LEF consist only of positive portfolio fractions. The solid square represents the Kelly 
portfolio as follows from Eq. ([9|) (again, short selling and borrowing are forbidden). We see 
that EF and LEF are close to each other and thus from the practical point of view they do 
not differ. 

Finally, let us discuss a useful simplification which allows us, in some cases, to reduce the 
time-consuming numerical computations. By differentiating Eq. ()16p we obtain the condition 
for the optimal portfolio fractions 

K^^-)"^'Kt)"^^ = " (.: = V^.iV). (IT) 

When the parameters of all assets fulfill the condition mi,Di <^ 1, we can use approximations 
for E{Ri/Wi) introduced in Appendix^ The first term can be evaluated more precisely using 
In Wi = ln(l + Rp) « Rp. Hence E{Ri In Wi/Wi) « E{Rj) + Y^j^i Qj E{Ri)E{Rj), where 
Ri = e''' — 1 is the return of asset i. Furthermore, for m^, Di <^ 1 we have E{Ri) mi + Di/2 
and E{Rf) ^ Di. As a result we obtain the equations 

2g,A + 2{mi + A/2) ^ qj{mj + Dj/2) + 71 



Di , 

"ii + ^ (1 - %) 



+ 72 = 0, (18) 
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Figure 7: Comparison of two efficient frontiers. The solid line is the classic EF given by 
Eq. ([5]). The dashed line is obtained by fixing E{lnWi) and by numerically minimising 
fluctuations of In W"i around this value. This we call LEF, its portion with positive qi values 
is highlighted in thick gray. The small filled circles represent the three individual assets that 
compose the system. The filled square is the Kelly-optimal portfolio (no borrowing, no short 
selling). 

where i = 1, . . . ,N and the values of 71 and 72 are fixed by the constraints Xli^i Qii^i + 
Di/2) = vp, YliLi Qi ~ 1- This set of nonlinear equations allows us to approximately find 
Logarithmic Efficient Frontier. In comparison with a straightforward numerical maximisation 
of Eq. (jl6p (involving numerical integration of £^(ln Wi) and £^((ln Wi)'^)), a substantial saving 
of computational costs is achieved. 

4 Concluding remarks 

In this work we investigated the Kelly optimisation strategy in the framework of a simple 
stochastic model for asset prices. We derived a highly accurate approximate analytical formula 
for the optimal portfolio fractions. We proved that in the limit of small returns and volatilities 
of the assets, the constrained Kelly-optimal portfolio lies on the Efficient Frontier. Based on 
the obtained analytical results, we proposed a simple algorithm for the construction of the 
optimal portfolio in the constrained case. We showed that since in the investigated case 
of lognormal returns, Kelly's approach forbids short positions and borrowing, only a part 
of the available assets is included in the optimal portfolio. In some cases the size of the 
optimal portfolio is much smaller than the number of available assets — we say that a portfolio 
condensation arises. In particular, when the distribution of the mean asset returns is wide, 
there is a high probability that only the most profitable asset is included in the Kelly-optimal 
portfolio. 

The Mean- Variance analysis is a well-established approach to the portfolio optimisation. 
We modified this method by replacing the averages E(Wi) and E{Wi) with the logarithm- 
related quantities E(lnWi) and ii^((ln 1^1)^). These are less affected by rare events and allow 
to capture the typical behaviour of the system. As a matter of fact, the difference between 
the traditional M-V approach and the modification proposed here is very small and does not 
justify the additional complexity thus induced. 
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A Main approximations 

Our aim is to approximate expressions of the type E{g{ri)), where r] follows a normal distri- 
bution f{r]) with the mean m and the variance D. For small values of D, this distribution 
is sharply peaked and an approximate solution can be found expanding g{rj) around this m. 
This expansion has the following effective form 

ef 1 1 

giv) = 9{m) + -{r] - mfg^'^'> (m) + ^ - mfg'''^^ (x) , (19) 

for some x G [m, rj]. Here we dropped the terms proportional to (rj—m)^ with an odd exponent 
k, for they vanish after the averaging. If we take only the first two terms into account, we 
obtain 

/oo 
9iv)f{v) dv ~ 9im) + - g^^Hm). (20) 
^oo ^ 

This approximation is valid when the following term of the Taylor series brings a negligible 
contribution A. We can estimate it in the following way (x E [jn-, ??]) 

A = y ^' 9^^Hx)f{r])dv< J ^' Mf{rj)d7j = —^. 

Here by M we label the maximum of [^^^^(t/)! in the region X where f{r]) differs from zero 
considerably, e.g. = [m — 2D,m + 2D]. Since g{x) has no singular points in a wide 
neighbourhood of m, its fourth derivative is a bounded and well-behaved function. Thus M 
is finite and A vanish when D is small. 

In particular, in this work we deal with functions of the form (^(T^j) = (e'' — l)/[l-|-g(e'' — 1)]. 
If we use Eq. ()20p with this g{i]), approximate 1 + q{e'^ — 1) in the resulting denominators by 
1, e"^ by 1, and e™' — 1 by m, we are left with 

E{g{rj)) ^m + D{l-2q)/2. (21) 

We widely use approximations of this kind to obtain the leading terms for the optimal portfolio 
fractions in this paper. 

B Procedure for correlated asset prices 

So far we have considered uncorrelated asset prices, undergoing the geometric Brownian 
motion of Eq. ([T]) . Obviously, this is an idealised model and real asset prices exhibit various 
kinds of correlations. In order to treat correlated prices we employ the covariance matrix S 
to characterise the second moment of the stochastic terms E{{r}i — mi){r]j — rrij)) = Sij. The 
uncorrelated case can be recovered with the substitution Sij = SijDi. 
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Again, we would like to find an approximation of the term E(^g{r])) = J g{r])f{r])dr]. 
Here f{r]) is the probability distribution of r/ and g{r]) is the function of interest. Notice that 
the correlations impose the use of vector forms for all the quantities of interest. The Taylor 
expansion of g{r]) around m, Eq. (I19p in the uncorrelated case, takes the form 

g{rj) = g{m) + Vg{m) ■ {rj — m) + -(r/ — m)^\/{m){rf — m) + . . . . 

Here V(m) is the matrix of second derivatives of the function g{rj)^ calculated at the point 
r] = m. Now we can proceed in the same way as before 

N 

E{giv)) ~ g{m) f{r])dT] + ^dig{m) {r]i-mi)f{r])dr] + 
•' i=i •' 

1 ^ f 

1 ^ 1 
= g{m) + -Y,S^,V,,=g{m) + -TT{Sy). 

In the last line we used the symmetry of S. For given g{ri), m and S, we can now solve the 
equation E{g{r})) = 0. In particular, these approximations can be cast into Eq. ([6]), which 
can then be treated as in the uncorrelated case. 
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